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I. INTRODUCTION 

The main cosmological constant problem [l|, |2| can be phrased as follows: why do the 
quantum fields of the vacuum state not naturally produce a large (positive or negative) value 
for the cosmological constant with an energy scale of the order of the known energy scales of 
elementary particle physics? 

An ideal solution would be to compensate dynamically any cosmological constant there 
may be. In equilibrium, such a compensation appears to be impossible with a constant 
(spacetime- independent) fundamental scalar field [1|. Partly for this reason, Dolgov |3|, 
^ has proposed using nonconstant higher-spin fields, notably a nonconstant vector field. 
He presented a remarkably simple cosmological model with a single massless vector field 
Aa{x), which allows for the compensation of a cosmological constant A of a particular sign 
with Minkowski spacetime appearing asymptotically as an attractor of the dynamical field 
equations. However, a serious flaw of this compensation-type solution to the cosmological 
constant problem was pointed out by Rubakov and Tinyakov [5|, namely, that the resulting 
Minkowski spacetime (with a vector-field background canceling A) implies an unacceptable 
modification of the standard Newtonian gravitational dynamics for small systems. 

In this article, we present a specific model with two massless vector fields, Aa{x) and 
Bp{x)^ which evades the above-mentioned flaw with the local Newtonian dynamics. Inspi- 
ration for this model was obtained from previous work by Volovik and one of the present 
authors on the so-called g-theory approach |6|-l9| to the main cosmological constant problem 



(a one-page review of g-theory can be found in Appendix A of Ref. 10| and a ultrabrief 
summary will be given in the Footnote of Sec. IIIAp . In Ref. [9|, in particular, it was real- 
ized that the Dolgov theory actually provides a generalization of g-theory, with the genuine 
g-theory appearing asymptotically. Therefore, the insights of g-theory can also be applied 
to Dolgov-type vector-field models and be used to overcome the Newtonian-dynamics flaw. 

II. MINKOWSKI ATTRACTOR FROM A VECTOR FIELD 
A. Generalized Dolgov model 

Our starting point is the vector-field model presented by Dolgov |3|, y] (related aether- 
type theories have been discussed by, for example, Jacobson ll|). Here, we extend the 
previous analysis of Ref. |9|, in order to compensate both positive and negative cosmological 
constants in a single model. 

The effective action of the massless vector field Aa{x) and the metric gapix) is taken to 



be the following {h = c = 1): 

S,s[A,g] = - jd'xV^(—^R[g]+e{Q[A,g])+AY (2.1a) 



Q[A,g] = ^A„;^A";^ (2.1b) 

where e is an appropriate function of the variable Q (semicolons in the definition of Q denote 
covariant differentiation), R is the Ricci scalar, Gn is Newton's gravitational constant, and 
A is the effective cosmological constant. The above action generalizes the one of Dolgov [^ 
which has e = —r]oQ^ for ?7o = ±1. As mentioned in Ref. j^, the consistency of such a 
massless vector field Aa at the quantum level (e.g., unitarity) needs to be investigated, but 
this issue lies outside the scope of the present paper which is primarily concerned with the 
classical dynamics of the metric and vector fields. Incidentally, this vector field A^ is not a 
gauge field, so also its masslessness needs to be explained, but, again, this issue will not be 
addressed here. 

It may be useful to present a simple example of a bounded function e{Q), which gives a 
unique equilibrium value Qq for each value of the cosmological constant A: 

= I ^-- (l - VI - (Q - Qm)V(AQ)2) , for Q G (Q,n - AQ, Q^ + AQ) , 
I , otherwise , 

for Qin > AQ > and a constant Cmax > 0. The corresponding gravitating vacuum energy 
density, given by e{Q) = t{Q) — Qdt{Q)/dQ according to Ref. [6|, descends monotonically 
from +00 to — oo as Q runs from Q-^ — AQ to Qm + AQ. This behavior oi'e{Q) indeed allows 
for the compensation of any value of A by a unique equilibrium value Qq [see, in particular, 
d2lbll below]. 

As suggested by Dolgov |3|, 1^, the following isotropic Ansatz can be taken for the vector 
field Aa{x) in a spatially flat Friedmann-Robert son- Walker (FRW) universe: 

Ao = Ao(t) = Vit) , Ai = ^2 = A3 = , (2.3a) 

(gap) = diag(l, -a{t), -a{t), -a{t)) , (2.3b) 

where t is the cosmic time and a{t) the FRW cosmic scale factor with Hubble parameter 
H = {da/dt)/a. 

The reduced field equations are given in App. |S1 With appropriate boundary conditions 
(consistent with the Friedmann equation), numerical solutions have been obtained. These 
numerical solutions show that, for either sign of the cosmological constant A, there exists a 
finite domain of boundary values V and dV/dt at t = tstart, which give the same asymptotic 



solution for t — )■ oo: 

y(t)~(go/2)t, H{t)^l/t. (2.4a) 

The particular value Qq entering the dynamical solution fl2.4ap precisely cancels the effects 
from the cosmological constant [g], 



A + 



'«)-«^ 



= A + ?(Qo) = 0, (2.4b) 



as the fundamental dynamic variable A"p = A^^.^ = V/3 A°' approaches the Lorentz-invariant 
tensor structure J9| characteristic of g-theory/ 



A'^.ix) 



_ 1 

equil 2 



Qo <5"/3 ■ (2.4c) 



This shows that Minkowski spacetime can appear asymptotically as an attractor of the 
dynamical equations considered, independent of the sign of the cosmological constant (figures 
similar to Fig. 2 of Ref. [9| have been obtained but will not be given here). 

B. Flawed Newtonian dynamics 



Rubakov and Tinyakov [5j considered the quadratic action of small changes in the fields 
away from the attractor solution (I2.4ap : 

A^{x) = Af{x) + v^{x) ~ (tgo/2) 5l + v^{x) , (2.5a) 

9ap{x) = Qa^ix) + hc^fsix) ~ r]ai3 + haf^ix) , (2.5b) 

where the perturbed fields are distinguished by a hat in order to avoid confusion later on. 
From (12. ip and (12. 3p . they obtain the following structure of the field equation for the metric 
perturbation ha/si^x): 

( ^ >, (GR) ^ 

(SvrG^)-! I" 92/1 "I +{A,f -a^/i" =re.t, (2.6) 



^ Very briefly, q-theory aims to give the proper macroscopic description of the Lorentz-invariant quantum 
vacuum where a (Planck-scale) cosmological constant A has been canceled dynamically by appropriate 
microscopic degrees of freedom (see the original article [6] for further details and also the brief review [10|). 
Typically, there are one or more of these Lorentz-invariant vacuum variables (denoted by g, with or without 
additional suffixes) to characterize the thermodynamics of this static physical system in equilibrium. A 
particular example of g-theory is given by the Lorentz-invariant vacuum variable Qo appearing in (I2.4cp 
with a value determined by the Gibbs-Duhem-type equilibrium condition (j2.4b[) . The issue discussed here 
is the dynamics^ namely, how the equilibrium state is approached. 



where the notation is symbohc with all spacetime indices omitted. The two occurrences of 
"9^ /i" in (12. 6p stand for different expressions, each involving two partial derivatives da, the 
Minkowski metric rjap., and the metric- field components hap- On the right-hand side of (12. 6p 
appears the energy-momentum tensor T^^^ of a local matter distribution. Note that (12. 6p 
for Aq = corresponds to the standard Einstein equation of general relativity (GR), which 
reproduces the Poisson equation of Newtonian gravity in the nonrelativistic limit. 

With Ao ~ ((5o/2)t, Qo ~ (SttGat)-! ~ (lO^^ GeV)2, and t ~ lO^^yr ~ (IQ-^s eV)-S 
the second term on the left-hand side of (12. 6p dominates the first term and ruins the stan- 
dard Newtonian behavior. This equation also suggests that the properties of gravitational 
waves are unusual compared to those from general relativity and, most likely, physically 
unacceptable |5[. 

III. MINKOWSKI ATTRACTOR FROM TWO VECTOR FIELDS 

A. Setup 

A possible cure for the flaw of Sec. Ill Bl uses two massless vector fields Aa{x) and Ba{x) 
with the following effective action: 



s-, 



eflf 



j d^x v^ Q (^Planck)' R + e(QA, Qb) + A j , (3.1a) 



Qa = ^jAa;p A--^P , Qb = ^jBa;p B-'^P , (3.1b) 

Epianck = (SttGjv)"^/^ ^ 2.44 X 10^® GcV . (3.1c) 

The Dolgov-type Ansatz for the vector fields Aa{x) and Bp{x) and for the metric Qapix) is: 

Ao = Aoit) = Vit) , A^=A2 = A, = Q, (3.2a) 

Bo = Boit) = Wit) , B, = B2 = B3 = 0, (3.2b) 

(fi'a/3) = diag(l, -a(t), -a(t), -a(t)) , (3.2c) 

where a{t) is again the cosmic scale factor of the spatially flat FRW universe considered. 
Solving the field equations from (13. lap for the Ansatz fields ( 13. 2p gives the explicit functions 
V{t),W{t), anda(t). 

For later use, we introduce dimensionless variables. Specifically, we replace the above di- 
mensional variables (and the variable X to be defined shortly) by the following dimensionless 



variables: 

{A,e,X,t,H} -> {X,e,x,r,h}, (3.3a) 

{Qa, Qb, V, W} ^ {qA, qs, V, w] , (3.3b) 

having used appropriate powers of the reduced Planck energy -Epianck without additional 
numerical factors. Moreover, |A| is considered to be of order unity. 

B. Main argument 

The action-density term e(Q^, Qb) for the two vector fields of the model (13. lap will be 
designed to cancel the effects of the cosmological constant A and to give a vanishing contri- 
bution to the field equation for the metric perturbation hap{x). Concretely, perturbations 
around the background solution from (13. lap and (13.20 give the following equation instead of 
(1^: 

( ^ >, (GR) ( ^ ^ ^ >, 

+ [^-^]as..p {^' "^'^" +^ "^^" + "^"} + [^+~^Ur., { "^"} = ^ext, (3.4a) 

[X-M = 0, (3.4b) 

L J asymp ' ^ ' 

\e-e] =0, (3.4c) 

[A + el =0. (3.4d) 

L J asymp ^ ' 

Equation (I3.4ap for the metric perturbation contains the asymptotic values of two basic 
quantities of g-theory p, |8[, namely, the inverse of the vacuum compressibihty (denoted by 
the Greek capital letter 'chi'), 

^-, ^ Q. d'<QA. Qb) , dMQA, Qb) 2 Q^ q^^MQa^QbI 

^^ dQAdQA ^"^ dQBdQs ^^ dQAdQs ^ ' 

and the thermodynamically active (and gravitating) vacuum energy density. 

Physically, conditions (]3.4bp and (I3.4cp can be interpreted as having a perfectly soft and 
flexible medium (isothermal compressibility X = —V~^ dV/dP = oo), which does not affect 
the metric perturbations. But, for the moment, we are only interested in finding a working 



model and follow Newton's advice, ''Hypotheses non fingo" |l2 |. 



The derivation of fl3.4ap proceeds in five steps. First, consider the second-order variation 
of the act ion- density term e{QA, Qb), 



Q'a 



(P 



+ 



(IQa dQA 
2 QaQb 






+ 



dh 



dQAdQs 



Qa 
6Q'l' 6Q 



Ql 

(1) 

B 



d\ 



5Q^^ 5Q%^ 



dQs dQB 
de 



Qa Qb 



+ 



Q. 



dQ, 



Qb Qb 

(2 
A 



SQ''' 



Q. 



+ 



Ql 



de 



dQ 



B 



SQl 
Q 



B 



(3.6) 



where all factors in square brackets are evaluated with the background solutions V{t), W{t), 
and a{t) from (13. lap and (13. 2p . 

Second, observe that all factors 5Qx /Qx in the above equation have the same structure 
and so do the factors SQx /Qx- In terms of the perturbative fields v, w, and h [the definitions 
of Va and hai3 have been given in (12. 5p . the one of Wa is similar], (13. 6p becomes in a symbolic 
notation: 



Ql 



d\ 



+ 



Q 



dQA dQA 
de 



fiv, hf + 



r>2 d-'e 
^^ dQs dQB 



f{w, hf + 



2 QaQb 



dh 



dQA dQB 



A 



dQ, 



i{v, h) + 



Q 



de 



B 



dQ 



B 



i{w, h) 



fiv,h)fiw,h) 
(3.7) 



where the explicit expressions for the linear function / and the quadratic function i can be 
obtained from the results given in App. [Bl 

Third, assume certain (anti-)symmetry properties of e{QA, Qb) and its derivatives and 
also the existence of a Dolgov-type asymptotic background solution [both assumptions are 
satisfied by the specific example of Sec. IIIICI for 6 = 0]. Then, it can be shown that the 
resulting equations for v and w have an identical solution, provided the matter perturbation 
is localized. The implication is that the first three terms in (13. 7p combine and so do the last 
two terms. Indeed, a direct calculation gives: 

d'^e „o d'^e „ „ „ rf^' 



Q\ 



+ 



Q 



dQA dQA 
de 



+ Ql 



dQB dQs 



+ 2 Qa Qb 



dQA dQB 



§2 dh dh + gi h dh + go h h 



+ Q 



de 



B 



k2 dh dh + ki h dh + ko h h 



dQA dQB _ 

= [X~^] lg2dhdh + gihdh + gohh\+[e-e] ^2 dh dh + kih dh + kohh\ , (3.8) 

with a symbolic notation for the prefactors gn and kn- 

Fourth, consider the further (standard) contribution to the quadratic action, which follows 
from the variation of the metric entering the spacetime measure of (13. lap . Specifically, this 
contribution is 



[A + e] {khh^ 



(3.9) 



again with a symbolic notation for the prefactor Iq. 

Fifth, make the necessary partial integrations in fl3.8p . while assuming X~^ and (e — e') 
to be constant on the macroscopic length scales considered, and add the contribution (13.91) . 
The resulting quadratic action then gives the linear field equation fl3.4ap . Note that, in the 
last term on the left-hand side of fl3.4ap , the asymptotic value of e has been replaced by the 
one of ?, in agreement with fl3.4cp . 

It now remains for us to present an Ansatz for e{QA, Qb) with appropriate symmetry 
properties and with both [X~^] and [e — e'] vanishing identically (that is, purely by algebra). 
This will be done in the next subsection. 

C. Specific model 

Following up on the general discussion of Sec. IIIIBt we now choose the action density 
e{QA, Qb) of ( 13. lap to be a particular rational function. Specifically, the dimensionless 
vacuum energy density e, the corresponding gravitating vacuum energy density e, and the 
corresponding inverse vacuum compressibility x~^ ^"^^ given by: 

(i^Pianck)^ ^ + {A.,, A-^P) [B^,, 5-;^) 6 + qWB 4 4 ^^ ' ^ ^ 

"?A aqB i^ + qAqB) ^^^B 

1 rf CL fZ- try CL (Z- (Jj C- 

X = qA-} — 1— + ?B3 j— + '2qAqB 



_4^ y----B^-jy.A .B> ^ _20 5 __ + 0(5^), (3.10c) 



A uqA dqs dqs dqA dqs 

{5qlq l-3 6) {q\-q%) _ ^^ ^ e 
{S + q\q%Y iWb 

for a positive infinitesimal 5 and dimensionless variables from (13. 3p . The last steps in the 
above three equations give the leading order in 5. It is certainly possible to set 5 immediately 
to zero in (I3.10p . but we prefer to keep 5 explicit in order to clarify two technical points 
later on, regarding asymptotes and stability. It needs, however, to be emphasized that the 
actual model function e{QAi Qb) for the action (13. lap is the one obtained from (I3.10ap with 
(5 = exactly, so that Eqs. (]3.4bp and (I3.4cp hold identically. 

From the Dolgov-type Ansatz (13. 2p . the following ordinary differential equations (ODEs) 



for v{t), w{t), and /i(r) are obtained (cf. App. |A|) : 

(v + 3hv-3h'^v] 

\ J QAdqA 

( , 2 '^ de 

( w + 3 hw — 3h w 

2h + 3h'^ = X + 



de , d / de 

dT\ QAdqA 

. d / de 
qsdqB drKqBdqs 



gyl=\/ir+3fi??J^, qs=\/«''^+3 h'^ w'^ 



, (3.11a) 
0, (3.11b) 



d / 2 de 

e{qA, Qb) --ry'^^ — 3" 
dr \ qA dqA 

d / r, de \ r, de 

dr \ qs dqs 



de 



+ v^ 



qAdq/ 



+ w^ 



qBdq_ 



B 



(3.11c) 



where an overdot stands for differentiation with respect to r. The corresponding Friedmann 
equation is given by 



3 /i^ = A + 



e(gA, qs) 



(3.12) 



qA=\JiP'^Z h'^ v'^ , qB=\/iJj'^+'i h'^ w'^ 

The boundary conditions for f (r), v{t), w{t), w{r), and /i(r) at r = Tstart niust satisfy (I3.12p 
with a nonnegative right-hand side. Observe also that, just as in (12.4b|) for the generahzed 
Dolgov model of Sec. Ill At the right-hand side of (I3.12p can be nullified by an asymptotic 
solution with the appropriate constant value for the ratio of the auxiliary variables qA and 

qB- 

The asymptotic behavior of the solutions of (13.1 ip is rather subtle. Mathematically, the 
order of limits 5 J, and r — t- oo is important. Physically, we take a fixed extremely small 
value of S and consider only "modest" values of the dimensionless cosmic time r: 



10 



-10^ 



r < 10 



60 



(3.13) 



where the last inequality includes cosmic times up to the present age of the Universe in units 
of ipianck = iZ-S'Pianck- It is, of coursc, possiblc to take a less radical value for 6, but the one 
in (I3.13P dispenses with some unnecessary discussion later on. 

For appropriate boundary conditions at r = Tstart = 0(1) and small but finite values of 
6, the solutions of (13. lip have the following asymptotic behavior for r — ?■ oo: 



V ^ kr^ , w ^ It^ , h ^ riT 



(3.14a) 



k''/f = v/(A/2)2 + l-A/2, (3.14b) 

= p{p-l)-3np + 3n^ , (3.14c) 

= S[p'^ -16np + 5n{4: + 3n)], (3.14d) 

where the parameter ratio (]3.14bp follows from (I3.12p . the relation (I3.14cp from (13.1 lap , 
and the relation (I3.14dp valid for p > 1 from the pressure terms in (13.1 Id) . Equations 



fl3.14cp and (13.14dp for (5 7^ give two sets of values (n, p) with p > 1. One set has values 
(n, p) ^ (0.6480, 2.424). But it is the other set, with values 



2 

n = 

183 



1 Q73771 

83 + V 14441 cos ( - arccos ■ 



3 (14441)3/2 



2.152, (3.15a) 



3rl+ 5 
p = An _ ?^ 3.655, (3.15b) 

13r2 — 1 

which will turn out to be relevant for the numerical results to be presented shortly. As 
mentioned in Sec. IIIIBt the Dolgov-type asymptotic solution (I3.14ap enters the derivation 
of (I3.4ap . as do the (anti-)symmetry properties of (I3.10al) and its derivatives for 5 = 0. 



D. Numerical solutions 

For A = ± 2, (5 = 10^^*^, and boundary conditions in an appropriate domain, the numerical 
solutions of the reduced fields equations display an attractor-type behavior (Fig. [T]) with 
V (X T^, w cc tP, and h r^ nr^^ for coefficients p and n from 03.15p . Different A values and 
different boundary conditions are seen to give an identical asymptote h ~ ?Tr~^ [remark 
that the normalization of v{t) is irrelevant for this asymptote; what matters is the constant 
ratio Qa/qb]- Within the numerical accuracy, the same results have been obtained for 5 = 0. 
The issue of the allowed boundary conditions is, however, more complicated and a complete 
analysis is left for a future publication. Another topic for future investigations is the possible 
cusp-like behavior at r ~ 1.6 suggested by two A = — 2 solutions hij) in Fig. [H 

The h panels in the third column of Fig. [1] show the main result: the approximately 
constant Hubble parameter /i of a de- Sitter-like universe at r ~ 1 changes to /i ~ r~^ for 
r ^ 1, so that a Minkowski spacetime (with /i = 0) is approached asymptotically. Moreover, 
the flat spacetime obtained for small but nonzero 5 has inverse vacuum compressibility 
Xasymp = as p > 1. The actual p value from (I3.15P even gives lim,-_j.oo t^X'^t) = 0, 
as required by (I3.4ap . Of course, [x"^(t)] vanishes identically for 5 = 0, and so does the 
quantity [e(r) — e(r)]. Finally, the total gravitating vacuum energy density [A + e(r)] is 
found to drop to zero as r~^, in agreement with (13.121) . As remarked already in Eq. (5.16) 
of Ref. |7| (and reiterated in the review [lQ[), the present value of this quantity [A + ?(tnow)] 
would be of the order of the experimental value A'^'^^p^ ~ (meV)^. 

E. Remarks 

Several points about the proposed A-cancellation mechanism of this section are to be 
noted: 
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FIG. 1: Numerical solutions of ODEs (13.1ip with dimensionless cosmological constant A = ±2, 
energy density function (|3.10ap . and model parameter 6 = lO"^'^. The following auxiliary functions 
are obtained from v{t), w{t), and h{T): qa = ^/ifl + 3 /i^ f ^ and qB = \/'uP' + 3/1^ w^. 
Top row (A = +2) : The boundary conditions are {w(l), w{\)\ = {8/(3/2 - 1/100 + r/25), 8} and 
{ij{l), w{l)} = {(3/4 + s/2)/(3/2 - 1/100 + r/25), (3/4 + s/2)} for integers r = ±1 and s = ±1. 
The corresponding values for /i(l) follow from (j3.12p . The dashed lines in the plots refer to the 
lowest value of v{l) coming from s = —1. The scaling of the v{t) plot uses the function (T-(-(r) = 
[1 + 35 (r - 1)2]/ [10 + 100 (r - 1)^ + (r - if+P] and the scaling of the /i(r) plot uses r/n with 
exact parameters {n,p) from (13.15p . 

Bottom row (A = —2): The boundary conditions are {v{l), w{l)} = {6, 6/(2 + r/5)} and 
{vil), w{l)} = {(1/2 + s), (1/2 + s)/(2 + r/5)} for r = ±1 and s = ±1. The scaling of the 
v{t) plot uses the function o-_(t) = [2 + 12 (r - l)2]/[30 + 120 (r - 1)^ + (r - 1)2+p] and the 
scaling of the /i(r) plot uses r/n with exact parameters {n,p) from ()3.15p . 



(i) With p > 1 from (I3.15p . a nonstandard form of g-theory is obtained asymptotically, 
having growing individual values ^aIt) ~ r^^^ and ^^(r) ~ r^^^ but a constant ratio 
Qa/qb- This behavior allows for both the cancellation of the cosmological constant A 
and having lim^^^oo X"^(^) = Xasymp = ^■ 

(ii) Even if p were equal to unity [which, with n = 1, is also a possible solution of the 
reduced field equations (13. lip ], the present values of X~^ and (e— ?) would be negligible 
for the S value displayed in (IS.lSp . bringing (I3.4ap extremely close to the standard 
Newtonian result. 
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(iii) An entirely open issue is the question of stability (cf. Ref. |6[), where f l3.10cp becomes 
asymptotically +20 6 A/(g^g^), which is only positive for the case of A > (5 being 
positive by definition). The possible instability of the A < solution may be consistent 
with the fact that the numerical A = —2 solution of Fig. [1] has been found to become 
ill-behaved for 6 ~ 10~^ (i.e., divergent at finite values of r), whereas the numerical 
A = 2 solution for 6 ~ 10~^ remains unchanged compared to the 6 = 10^^° case. 

(iv) In the very early universe, i.e., far away from the asymptote, the perturbation equation 
fl3.4ap differs from the standard Einstein expression. This different equation may lead 
to new effects for the creation and propagation of gravitational waves in the very early 
universe (assuming the model of this section to be physically relevant). The main focus 
of the present article is, however, on the Newtonian physics in the final equilibrium 
state of the universe. 

IV. CONCLUSION 

The fundamental question addressed in this article is whether or not a vector-field 
model [sl, |j| allows for the dynamic cancellation of an arbitrary cosmological constant A 
without spoiling the local Newtonian gravitational dynamics |5|. The answer found is af- 
firmative, even though the final FRW-type universe obtained {H ~ 2t^^) does not quite 
resemble the actual Universe of our recent past {H ~ nt^^ for n changing from 1/2 to 
2/3). The important point is that, as a matter of principle, it is possible to evolve from an 
initial de- Sitter- type universe [with a cosmological constant |A| ~ (-Epianck)^] to an asymp- 
totic Minkowski spacetime [with A^fj = A + ^{Qaq/Qbo) = and standard local Newtonian 
gravitational dynamics]. 

It is clear that the explicit vector-field example of Sec. IIIICI can be generalized. It may 
even be possible to appeal to higher-spin fields, perhaps the well-known threeform gauge 
field (cf. Refs. [g, |9| and further references therein). The most important task, however, is 
to establish the consistency of this type of vector-field model and to discover the underlying 
physics. 
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NOTE ADDED 

The present article considers a particular A-cancellation vector-field model which, asymptot- 
ically, has the standard Newtonian dynamics on small scales but not an acceptable Hubble 
expansion on cosmological scales, as noted in the first paragraph of Sec. IIVI This article is, 
in fact, the first of a t rilog y of articles. 

The second article 13|] of the trilogy considers a different model which, asymptotically. 



gives the standard radiation-dominated FRW universe with H = {1/2) t ^ but, most likely. 



not the standard loca! 



N 



ewtonian dynamics. 

The third article IJ] of the trilogy considers a final model (combining the essential 
features of the two previous models) which, asymptotically, has both the standard local 
Newtonian dynamics and the standard radiation-dominated FRW universe. This last arti- 
cle also gives a mathematical discussion of the attractor-type behavior found in the three 
different vector-field models considered. 

A further article [l5|, a direct follow-up of the present one, discusses the possibility of 
having an early-universe phase with inflation and a late-universe phase with a dynamically 
canceled cosmological constant A. 

Appendix A: Field equations 

The action (12. lap gives the following field equation for the vector field Aa{x): 

V°(CV„A;3)=0, (Al) 

in terms of the function ({Q) = e'{Q)/{2Q), where the prime denotes differentiation with 
respect to Q. For a spatially fiat FRW universe, ( lAll) reduces to 



C [a"a„ + 3Hdo- 3H^ + C"^C'"5a] Ao-[2CHd^ +H C'^] A^ = 0, (A2a) 



C 



d'^da + Hdo-H- 3H' -C CH + C~ C'^^a Aj + [2CHdj + H Q] Aq = 0, (A2b) 



where, in this appendix, an overdot stands for differentiation with respect to the cosmic time 
t and H is the Hubble parameter defined as a/a. Furthermore, the quantity C'" denotes d°'(, 
the index a runs from to 3, and the index j runs from 1 to 3. 

The energy-momentum tensor T^piA) is obtained by varying the action (12. lap with re- 
spect to the metric ga/3 : 

T„^(A) = T^„(A) = e(Q) ga^~2C [A,.^ A'^ + A^,^ A^] 

+V^ [C {AaA^;p + Ap A^,^ + A^ Afi,^ + Afi A^,^ - A^ A^,^ - A^ Ap,^)] , (A3) 
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where Aa-^ denotes the covariant derivative V^Aa and similarly for other tensors. An alter- 
native form of this energy-momentum tensor is 

T^fsiA) = [eiQ) -CQ'] gap - 2CTX''^^"^(^) 

+ (V^C) [Aa A;/3 + ^/3 A;" + ^» Mj + ^/3 Aan - A (^«;/3 + Ma)] , (A4a) 

-'■a/S \^) ~ ~'^Q 9al3 + ^a;7 A'^ + A^-a A.^ 

--V^[A^A^,p + Ap A^.^ + A^ Ap,^ + Ap A^,^ - A^ A^.p -A^Ap,^], ( A4b) 

where TT^ ratic^^^-j ^^^^^^ with expression (7) of Ref. [^ for r/o = +1. 
The isotropic Ansatz (12. 3p reduces (]A2aP and (lA2bp to a single ODE, 

io + fsif + C/C) ^ - 3if 'Ao = 0, (A5) 



assuming C, to be nonzero. Note that C, in the above equation is a function of Aq. The 
implication is that flASP is, in general, nonlinear in Aq. 

Similarly, we can find the Ansatz energy density p{A) [from the definition Tq^{A) = p{A)] 
and the isotropic pressure P{A) [from the definition T- '(A) = —P{A) 5*]: 

p{A) = e{Q) - Q ^ , (A6a) 

with 

Q2 = (io)2 + 3/7^2. (A6c) 

Finally, the isotropic Ansatz ( 12. 3 p reduces the Einstein field equations to the following FRW 
equations: 

3i7^ = 87rGjv[A + p(A)], (A7a) 

2iJ + 3i72 = 87rG^[A-P(A)], (A7b) 

in terms of the vector- field energy density and pressure from (1A6P . 

Appendix B: Quadratic perturbations 



Following the discussion of Ref. pj, we consider matter perturbations with timescales 
and lengths very much smaller than the cosmological timescale H^^ ~ 10^° yr and size 
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c/Hq ~ 10^^ m, defined in terms of tlie measured Hubble constant Hq ~ 75 kms~^Mpc~^. 
Tfiese matter perturbations are considered to be relevant to the local Newtonian dynamics. 
Perturbing around the Dolgov-type solution f l3.2p . the second-order variation of the La- 
grange density (13. lap of the two vector fields reads: 



£(2) = C 



with 



£ 



(2) 



(2) 






-B 



2Q, 



AB ' 



1 rfe 



dQA \QAdQA 



\ ^a,P^r,& + 



de 
dQA 



r'9^' 



X {bA^,fi bA^,, - 2 bA^,fi bY% Ao + bY% bY% Ag) 



(Bl) 



(B2a) 



£ 



(2) 
B 



2Q 



B 



d 



1 de 



dQB xQBdQs 



Qa;l3Qr,& + 



de 



X ( 5Ba-p 6B^.s 



2 5B^,p5T%B^ 



dQ 



g^lgPS 



B 



(5r°^ (5r°5 Bl 1 , 



(B2b) 



£ 



(2) 
AB 



QaQi 



£e 



dQAdQs 



j^a;Pj^-y;5 



X (^5A^,^ 6B^,s - 6A^,p 5r°^ Bo - 5B,,s SVl^ A^ + 51°^ dV^, A^ B^) , (B2c) 

where 5Aa{x) and 5Ba{x) are the vector perturbations and (5r°^(a;) = (1/2) [/?.oa,/3(a^) + 
hoi3,a{x) — hai3fl{x)] coutaius the metric perturbation hai3{x). 

Quadratic terms of order H hdh and H^ h"^ have been calculated but are not given 
explicitly in (]B2p . because they are subleading compared to the {d h)"^ terms shown [the 
relevant timescales for the Newtonian dynamics (e.g., in the solar system) are very much 
smaller than the cosmological timescale H~^]. Remark, finally, that the perturbation fields 
SAa, SBa, and hap entering (IB2p are denoted Va, Wa, and haj3 in Sec. IIIIBj see also the 
earlier definitions in (12.50 . 
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